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ABSTRACT 


Damping is one of the most crucial and complex parameters 
that determine the performance of a rotor. Rotor damping is a 
result of internal < intermolecular friction i.e. material 
hysteresis) and external (Coulomb friction at supports and 
gas/air environment resistance) energy dissipation mechanisms. It 
is essential to make a distinction between contributions from 
these two energy dissipation mechanisms since they affect the 
rotor behaviour in different ways. 

In this study an experiments. technique is developed for 
'Separate' estimation of the internal and external damping 
parameters of a rotor. 

A rotor test rig is fabricated. The transient polar plot of 
the system is experimentally obtained. Subsequently the known 
analytical equations governing the rotor whirl are modelled to 
make a distinction between the contributions from internal and 
external damping mechanism. 

Damping is analysed as a function of whirling mode, rotor 
speed and whirling amplitudes. 



NOMENCLATURE 


c 

e 

C. 

1 


X 5 y , X 5 y 




r 

Q 



Y 

o 


External damping coefficient 
Internal damping coefficient 

Displacement in vertical and horizontal directions 

Velocity is vertical and horizontal directions 

Initial conditions for displacement and velocity in 

vertical and horizontal directions 

Amplitude of whirl 

Amplitude of whirl at time t = O 

Cg^ /x ) = Non-dimensional error in vertical direction 
o 

<Q^/x ) = Non-dimensional error in horizontal directio 
^2 o 

(y /x ) = Non dimensional initial displacement in 
^o o 

horizontal direction. 


(x /x ) = Non dimensional initial velocity in vertical 
Q o ^ 

direction 


(y /x ) = Non dimensional initial velocity in 
horizontal direction 


X 


Y 


(x/x ) Non dimensional displacement in vertical 
o 

direction after time t 

(y/x^) Non dimensional displacement in horizontal 
direction after time t 



T 

P or 
a 



P 


OlV 


CA , . 

1 J 


(p X t) = Non dimensional time 

The damped natural frequency in rad/sec 

(6j/p) = Frequency ratio 

the rotational speed in rad/sec 

Influence Coef f icient . 
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CHAPTER I 


INTRODUCTION 


The need far higher speed, yet reliable operation of rotating 
iuachinery continues to increase. The key factor in development of 
high performance rotor is the ability to accurately predict the 
dynamic response and stability of the system in the design stage. 
Damping is one of the most crucial parameters that define the 
dynamic behaviour of the machine as well as its susceptibility to 
self-e:<cited vibrations leading to instability. 

The three primary damping mechanisms influencing the rotor 
performance are i) Viscous damping due to the surrounding gas 
medium, ii) friction damping at various component functions along 
the rotor, most frequently at the bearings and between the hub of 
the disc and the shaft, iii) material hysteresis or hysteretic 
damping which is caused by the intermolecular forces in the 
material arising as a result of the rate of deformation in the 
elastic shaft. Mechanism (i) is an external source of damping, 
while mechanisms <ii) and <iii) are categorised as internal 
sources. Unlike stationary machines where the total damping in 
the system is taken as the direct summation of the contributions 
from various damping mechanisms, in rotating machines these 


mechanisms affect the rotor behaviour 


in 


different ways. 
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Different sets of values of the internal and external damping 
parameters lead to either stable or unstable rotor behaviour, the 
unstable speed regime being defined by their relative magnitudes. 
Therefore, while in stationary machines it is sufficient to 
estimate the total effect of these damping mechanisms, in rotating 
machines it is essential that their contributions are known 
separately. 

Conventional methods of obtaining the damping values from the 
logarithmic decrement of free vibration signals of machines can 
not distinguish between the contributions from internal and 
external sources. In this study it is proposed to develop an 
experimental method for separate determination of external and 
internal damping values for a rotor. It is also proposed to 
define damping as a function of whirling mode, rotor speed and 
whirling amplitudes to refine the existing analytical models for 
stability predictions. 

1.1 LITERATURE SURVEY 

A number of authors have dealt the problem of unstable shaft 
behaviours due to the self excited mechanisms of internal damping 
at rotational speeds above the first critical. In these studies 
analytical models have been developed to study internal damping 
arising due to both material hysteresis as well as rubbing 
friction at various functions across the rotor. The stability 
criteria have been established by comparing the internal (i.e. 
destabilising) damping forces with the external (i.e. stabilising) 
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damping forces. Analytical expressions defining unstable rotor 
speed regimes in terms of the relative magnitudes of internal and 
external damping parameters have been obtained. 

The problem of unstable shaft behaviours due to self excited 
mechanisms of internal damping at rotational speeds above the 
first critical speed was detected as early as 1924 by Newkirk CIU 
and Kimball C2,31. 

Since then attempts have been made to study this tendency 
towards self excited vibrations. Ehrich C4,5] analysed the 
important role of dry friction and internal damping in inducing 
the instability in rotating systems. The rotational speed at 
which instability occurs is governed by the ratio of external 
friction to internal friction. 

Begg C63 developed an analysis showing that the regime of 
stability does exist and the induced whirling may propagate, 
persist or decay when the disturbance causing the whirl is 
removed, depending upon the frequency to damping ratio and upon 
the magnitude of initial disturbance. 

Bansal and Kirk C73 described an analytical method based on 
transfer matrix approach to calculate the damped critical speeds 
and instability threshold speed of multimass rotor bearing 
systems. The effect of bearing as well as bearing support 
flexibility on the stability of the rotor system were discussed. 

Zorzi and Nelson C8] through finite element approach 
investigated the effect of axial torque on the lateral dynamics of 
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rotor bearing systems- Glasgow and Nelson [91 presented the 
stability analysis of rotor bearing system and developed a method 
for significant reduction in the size of the overall system 
problems while retaining the essential dynamics characterist ics. 

1 . 2 PRESENT STUDY 

The analytical tools for prediction of rotor behaviour, 
discussed in the previous section, do however, require numerical 
estimates of the internal and external damping properties of the 
specific rotor- In the present study an attempt has been made in 
the direction of estimation of these damping parameters. 

The problem is defined and the derivation of the governing 
equation is briefly reviewed in Chapter 2- The free vibration 
solution to the equation is obtained and reformatted to describe 
the technique to extract information about the internal and 
external damping - 

Chapter 3 describes the experimental rig and subsequent 
investigations. The technique is validated and results are 
discussed . 

Chapter 4 summarises the present investigations and outlines 
the scope and possibilities of future attempts. 



CHAPTER II 


EQUATION OF MOTION AND RESPONSE 


In this chapter the derivation of the equation of motion for 
a rotor with a massless flexible shaft supported in bearings at 
its two ends and carrying a rigid disc at its midspan is reviewed. 
Such a model is popularly known as the 'Jeffeot rotor'. It 
represents the modal system obtained for any multimass rotor after 
modal analysis. Internal and external damping are both present in 
the model. The solution to the equation of motion is presented. 
The solution is reformatted and the analytical basis for 
extraction of information about internal and external damping 
parameters is developed. 

2.1 EQUATION OF MOTION 

The Jeffeot model for the i-th mode of a rotor is shown in 
Fig. 2-1. Let x-y— z and be the stationary and rotating 
coordinate systems respectively <Fig. 2.2). The geometric centre 
of the disc is at c and the centre of mass is located at point G. 
For a rotational speed w the free-body diagram of the disc is 
described in Fig. 2.3. 


The equations of motion 

rax + (k? + h? ) cos cot - 

2 


= ma 


cos 


can be written as CIO, 113 

ikf) + hr)) sin 6>t + cx 

'jit - ma„ o" 

2 


sin at 


( 2 . 1 ) 


w 



Fig. 2.1 The Jeffcot Model for the Mode of a Rotor 
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my + (kf + h? ) sin o>t + 


= mSj to sin 


(kY) + hTQ) cos wt + cy 

U't - cos (j^t 


where , 


a^ = a cos ^ = a sin B 

1 o’ 2 o 


Defining radial vectors r 
r = C equations (2-1) 

r' -4- 2C r 20 ( r — 

e i ' 

In the above 



= K + ly, C = € + 
and (2-2) can be 

2 2 
ii*:^r> p*“r = a^ii-'*" 



k 


m 


IT} and 
combined 


e 


+ 


1 


( 2 . 2 ) 

( 2 . 3 ) 

noting that 
into 



i6u't 


e 


+^g 


(2*^> 


< 2 . 5 ) 


2.2 SOLUTION OF EOM 

The free vibration response of the above equation can be 
taken as 

r = A eK'p (iXt) (2.6> 

Substitution of in (2.4) gives 

- 2iX (C + C.) - + 2i6Y C. = 0 (2.7) 

e 1 1 

yielding , 

^ = i (C + C. ) ± C - (C + C. )^ - 2iu C. 

1,2 ei ei 1 

2 2 2 

Defining the damped natural frequency, = P - (C^ + C^), 

one gets 

X. ^ = i (C + C.) ± C P^ - 2iaj C. 

1,2 e 1 d 1 

iaiC. „ 1 /r) 

= i (C + C.) ± E ( P^ — 3 

e 1 d P^ P , 

d a 


9 


ojC 

i 2 

Since C. is a small quantity ignoring the term ( ) , one gets 


X 


1,: 


_ ^ ■ 

± P . + i (C + C . + - p ^ ) 

d e 1 


( 2 . 8 ) 


Hence the solution becomes; 


r = exp 


+ R 


[-< 
2 .XP [ - 


uC. 

^i - P 


. “I 

■^"1 


d 

<JC 


exp ( + ^ 


(C + C 
e 


. -1 

. + ) t exp ( - iP^t ) (2.9) 


For stability of the response (2.9) it is essential that 


C + C. 
e 1 


6jC. 

1 


( 2 . 10 ) 


It can be seen that the rotor will be stable at 

C 


6a < P . ( 1 

d 


C. 

1 


) 


( 2 . 11 ) 


2.3 EXTRACTION OF DAMPING PARAMETERS CC AND C, 5 FROM THE RESPONSE 

& i 

The technique for obtaining the damping parameters and 
from the free vibration response is formulated below. Rewriting 
the response of eqn. (2.9) 

oJC. 


= Rj exp - 
R^ exp 1^ - 


(C + C. 
e 1 


d 
<uC. 

(C + c. + 

e 1 P 


— ) t 1 exp ( +iP^t ) 

fi 


i- ) t j exp { -iP^t ) 


letting 


6jC. 


A = (C + C. ) , 

e 1 P . ’ 

d 


ioC. 

B = (C + C . + ) 

® ^ ^d 


and since 
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P=zx + iy, + i 

equation (2.9) can be written as 

—At 

r = (X + lY^ ) e (cos P.t + i sin P.t) 

11 a d 

*4“ (X„ 4- iY^) e (cos P.t + i sin P.t) (2-12) 

2 2 d d 

Separating the real and imaginary parts, we get 

X = (X, e“®^) cos P.t + {-Y, + Y.., ) sin P.t 

12 d 1 -i d 

(2.13) 

y = (Y. .e e ) cos P ,t + (X., -e - X., e sin P t 

■^12 d 1 d 

(2.14) 

In the above it should be noted that X^, X^, Y^^ and Y^ are 

constants depending on initial conditions. 

Differentiating (2.13) and (2.14) with respect to t, 

• —At — Rt 

x = (X, .e + X^.e ) (-P_, sin P^t) 

1 2 d d 

+ cos P^t {-AXj.e“^^ - BX.,.e“®^) 

+ (-Yj.e”^^ + Y.,.e~®^) (P^ cos P^t) 

+ sin P.t (AY..e~^* - BY^.e"®*) (2.15) 

d 1 2 

• —At — Rt 

y = CY, .e -t Y^-e > (~P^ sin P^t) 

■'I 2 d d 

-At -Bt 

+ cos P^t (-AY^.e - BY^-e ) 

d 1 2 

+ (X, .e~^^ - X^.e“®*) (P^ cos P.t) 

1 2 d d 

+ sin P^t (-AX^.e“^^ + BX^.e"®^) 


( 2 . 16 ) 
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Let the initial conditions be 


K = y = y^; 


and y = y at t = 0. 
o o 


Substituting in equations (2.13) - (2.16) and solving for 


and Y^, we obtain 


^ = 


X., = 


-((B-A) ;; + (B'-AB + 2P_.^) x + (A+B) y + 2P_. y > 

o do d’^o d o 


= 



t(B-A)^ + 4 P^^ 
d 

3 


C-(B-A) X 

o 

- 2P_, y + (A^ - AB 
d ■'o 

+ 2P ^) 
d 

X - P. (A+B) y > 
o d -^o 


C (B-A)^ + 4P_.^3 
a 



€2i< P. ■ 
a d 

- y (B-A) + X P . (A+B) - y 

^o o d -'o 

(B^-AB+2P ) > 

d 


C(B-A)^ + 4P_,^3 
d 



r o w P 

(. X r . — 

o d 

y (B-A) + P . X (A+B) + y 

•'q do ■'o 

(A^-AB+2P_,^)} 

d 


d 


Equations (2.13) and (2.14) can now be written as 


K = 1 , 5 - f C{(B-A) k + 2P_, y + (B^-AB + 2P_^) x 

[(B-A)^ 4P,-3 L o d ='o do 

d 

-At . . 

4- (A-4B) y >.e + C-CB-^A) x - 2P^ y 

d o o d o 

+ (A^ - AB + 2P_,^) X - P_ (A+B) y >.e“®^3 cos P_.t 

d o d ■'a d 

+ LC 2P_,.x - (B-A) y + P^ (A+B) x 

do •' o d o 

- (B^-AB+P.^) y >.e'’^^ + C 2P_, x -(B-A> y + P. (A+B) x 
d ■'o do “'o d o 

"7 2 — Ri- 

+ (A"'-AB+2P . ) y >.e 3. sin P.t 

d o d 


(2.17) 



C(B-A)^ + 3 


[e< - 2P^ K 


+ (B-A> y - P_,(A+B) X 
c ado 


■7 '? —At . . 

+ (B‘'-AB+2P_.^) y > .e + C 2P_, x -(B-A) y + P .(A+B) x 
d -^o do -^o d t 

+ (A^-AB+2P^^) y > e~®*3 cos P^-t + C-C(B-A) x + 2P_, y 

d d a d a 

+ <B^-AB+2P_.^> X + P_, (A+B) y > + tCB-A) x 

d o d o o 


* 9 ^ -Rt 

2P^ y -(A -AB+2P^^> x + P^(A+B) y > e 3 sin PA 
d o doo o 


d 

mM 

(2.18) 


Expressions (2.17) and (2.18) form a set of nonlinear 
algebraic simultaneous equations. It can be readily seen that 
knowing 

- the initial conditions x , x , y , y 

o D o o 

- the damped natural frequency pd 

- the response x, y after time t, 

the equations can be solved for A and B and hence the damping 

parameters C and C. . 

e 1 

Iterative Newton— Raphson technique in conjunction with 

experimental investigations can be employed to solve the equations 

and obtain C and C. . 

e 1 

Starting with assumed values of A and B the resultant 
response after time t can be calculated from equations (2.17) and 
(2.18) . 

Let the response thus calculated be denoted by x and y. 
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Let the actual response of the system after time t, as 
ibtained from experiment be x and y. 

The errors are 



The errors are minimised using Newton— Raphson technique and 
corrections • dA and dB are sought for the assumptions A and B. 
These corrections are obtained as 



Taking the corrected values of A and B the new assumptions the 

process can be repeated until the errors fall within specified 

<99, <9g dg Og 

limits. The detailed expressions for g, , g_, -xrr — , -s^ -. -et— , — 

1 2 oA aB OH oB 

used in the computational process are given in non-dimensional 


form in Appendix 1. 



CHAPTER III 


EXPERIMENTAL INVESTIGATIONS 


Ll THE RIG 

A portable rotor rig was fabricated and set up in the 
aboratory to carry out investigations. The rig is modular in 
Issign and can be employed for a variety of experiments on rotors. 

The rig consists of a shaft which supported in two identical 
all bearings. The bearing housings are mounted in the slots 
rovided in a base frame. Provisions are made in the base frame 
o carry shafts of different lengths- The shaft diameter can also 
e varied and correspond ing set of bearings can be supported in 
he same bearing housings- A number of discs can be mounted on 
he shaft as desired. The shaft is driven by a 50 W AC/DC motor 
hrough a flexible coupling. The motor is also mounted on the 
rame- A variac was used to control the input supply voltage to 
he motor and hence the rotational speed can be maintained as 
esi red - 

To set the rotor into transient whirling motion at any 
esired rotational speed, the system was inflicted with an impact 
oad. The transient motion can be readily recorded as a decaying 


lot . 
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2 INSTRUMENTATION 

The instrumentation Fig. 3.1 consists of two acce 1 erometers 
ounted on one of the bearing housing to pick up vibrations in 
he horizontal and vertical directions. The signals from the 
cceierometers are magnified in charge amplifiers and passed 
hroLigh filters (to filter the signal component with any required 
requency) and then displayed on a CRO- The CRO is connected to a 
rioter. In addition an FFT analyser is used for ident if ication 
f natural frequencies. 

The details of rig and instruments are given in Appendix 2. 
wo rotor conf igurat ions are considered 

i) Shaft with a single disc at its midspan (Fig. 3.2). 

ii) Shaft with two discs with unequal masses along its length 
(Fig. 3.3). 

.3 NATURAL FREQUENCY DETERMINATION 

Prior to experimentation for damping determination, the 
atural frequencies of the system need to be identified. This is 
ssential in order to filter out the response corresponding to a 
articular natural mode from the total response of the system as 
icked up by the accelerometers- 

The natural frequencies of the system are experimentally 
btained by subjecting the system to an impact load using an 
mpact hammer, so that it is set into free vibration. The free 

ibration response is displayed on the amplitude vs. frequency 
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6‘0 



Fig. 3.4 The FFT Plot for Single Mass Rotor 





Fig- 3.5 The FFT Plot for Two Mass Rotor 
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frame of the FFT analyser and the natural frequencies are 
ident i f i ed . 

Fig. 3.4 shows the FFT plot obtained for a single disc 

rotor. The natural frequency can be seen to be 44. S Hz. This 

value agrees well with the theoretical result obtained as follows. 

Idealising the rotor a single degree freedom system the 

natural frequency is 



For the configuration under consideration 

4SEI 2 

= 39584.067 kg /sec 
shaft ^3 ^ 

and m^. = 0.515 Kgs. 

disc 

yielding p = 44.2 Hz. 

Fig. 3.5 shows the FFT plot for the rotor with two discs. 
The natural frequencies can be seen to be 37.2 Hz and 128 Hz. 

The natural frequencies can be computed theoretically very 
readily by influence coefficients method for a two degree freedom 
system. 

Referring to Fig. 3.6 the influence coefficient matrix is 




ca) 


sia^^c a^/cc'S/^y^ c/ue ^cmS'S^ 



Ci>} IyG^‘^C'^^'cS>/y c/ure" ^ ^ocic/ /tj ciCJ^ S/ 


Influence Coefficients for Static Deflection of Shaft 



21 


where 


a. 


'21 


= 


-(8/81) L'^ 
6 El 

-(7/81) 

6 El 


The natural frequencies are given by C123 



a 


11 


m. 


+ ot_ 


m. 


± C (o( 


11 






ffl^) 


4m 


2 ^ 1/2 

"*12 


^ 1^2 


(ot 


11 '^22 



The theoretical natural frequencies are thus obtained as 32.5 Hz 
and 125 Hz and they can be seen to agree well with e;<periment. 

3.4 BALANCING 

The rotors need to be balanced before carrying out the 
investigations. Rigid rotor balancing was done for both the 
rotor configurations. Nominal amount of balancing was required. 
The details of balancing are given in Appendix - 3. 


3.5 DAMPING DETERMINATION 

The technique developed in Chapter 2 requires that starting 
with a set of initial values y^, x^ and y^ the modal free 
vibration response x»y of the system be known after time t. 

The rotor is run at the desired angular frequency iw by 
adjusting the motor input voltage. The rotor is then set into 
free vibrations by inflicting upon it a shock through the impact 
hammer. The filters are set to the desired natural frequencies so 
that out of the total system response picked up by the 
accelerometers only the component corresponding to set natural 
mode is displayed on the CRO. ^ - 
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The model response recorded simultaneously in the horizontal 
ind vertical directions at a rotor speed of 1200 rpm. are shown 
.n FiQ. 3.7 to These figures show the time response in the x 
[vertical) and y (directions) while Fig. 3.10 is a polar plot 
. .e. •! vs. y of the same signal. 

It IS known that for rotating components damping is not 
lonstant but a function of the mode of vibration, speed of 
'otation and amplitude of vibration C133. For a particular mode 
jf vibration and speed of rotation, damping will be a function of 
the amplitude of vibration. 

This phenomenon is manifested in the present study as well, 
IS will be seen later from the results. 

The iterative damping estimation procedure described in 
Chapter 2 is carried out for every cycle of decaying vibration 
signal. The damping being nonlinear in nature, for valid 
employment of the equation of motion described in the earlier 
ihapter, a practical approximation is made - that for a particular 
rotor speed and a specific mode of vibration the damping 
parameters remains practically constant over a cycle of vibration. 
Fhe damping parameters thus determined are construed to be valid 
only for the average polar amplitude of vibration during the cycle 
Linder consideration. Such approximations have been made in 
studies concerning rotating turbine blades 114,151 and is 
considered a pragmatic way of incorporating the dependence of 
damping on the amplitude of vibration. Hence assuming the damping 
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parameta-*'-3 .i::; a.’-fr i le the estimation process is 

repeated for ever'/ cvrle ai ".he decaying signal. The damping 
parameters ire then ’iir’.ctie * of the average polar amplitude 
during a c>cle of »ibratirjn. 

The d.i!Tiping estiin.it ion procedure is detailed m Table 1, for 
a typical signil (F'lg.Z.I! * ,at rpm, for the single mass rotor. 

Damping is evaluated bet'oesn joints 1 and 2, and 1' ^^d 2'. 

3.5 REStfLTS AND DI SCU^MT ON 

Computations as detailed in Table 1 '*jere repeated for five 

consecutive pair of amplitude in the decay signal- The plots of 

Amplitude ratio. F. tternal damping and Amplitude ratio internal 
damping for vinous rot.atianal speeds and modes for the case of 
rotor with signal mass at the midspan (as shown in Figs. o.l2 to 
3.16) and for rotor with two mass equally spaced between the 
bearings (as shown in Figs, 3.17 to 3.26 were made. 

It can be observed from the results that damping is dependent 
on the amplitude of vibration, rotor speed and the mode of 
vibration . 

For both the rotors it can be seen that the internal damping 
which is attributed to shaft material hysteresis increases 

with amplitude of vibration. This is in concurrence with the 

knowledge that material damping in metals increases with the 

amplitude of vibration. 




Ftg.3.15 Vorlaflon of dompfng wHh whlri ampITtude 
af 800 rpm for sTnglo mass rotor 


whirl cimpittwd# ratio ( r/ro) 


FT®. 3-1 4 V<irialiofi of damping with whirl ampll 
of 1200 rpm for singlo mass rotor 


whiri amplttu<l* ratio ( r/ro) 

H^molno with whlH omplltuds 


n9.3.15 Variation of damplno .T *1 1 rSnr 

ot 1600 rpm for slngis mass rotor 
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^ due to friction at supports and 

.iiTodyn im:r rr,;,*; ’ to decrease, though not 

i qn 1 f i u jn t i V witJ i i ? ■ ;ng amplitudes. More investigations 

under brlt-.:»r .. -;.ni », r? ’ » . v i ronment will be required prior to 

mal inq -iry on. hi .ton this phenomenon. 

Fo'‘ th ? two till* pa observed that both the 

et-ternal -ind ; r- 1 *» m t ! i'hpmg values are higher in the second 
mode. This i ^ m rsn 'f or n j with commonly observed phenomenon 

that damping i , the m the fundamental mode of vibration. 



CHAPTER 4 . 


CONCLUSIONS AND SCOPE FOR FUTURE WORK 


This technique for distinguishing between external and 

internal contributions to damping in rotors has been formulated 

established experimental ly- 

life rotors with 

This technique, on application to 

a. e-cr, ho oNtnected to fill a crucial gap of 

necessary modifications, can be exp 

. a. • 4- and 0 X t s TH a 1 damplnQ 

extractinQ information about in e 

^^chanisms separately in rotors. 

□nee information about internal damping characteristics of a 
newly manufactured rotor can be known by such a 
expected to be an important first step m design 
bearings whereby the bearing damping characteristics 

external damping for the rotor) can be appropriately modelle 
better control over rotor dynamics. 

Tbs laburstory big snd instrumsntstibn bomsver needs 
impruvement. The present study »es constreined by equipment 
availability. Instead of sensing the vibrations by bearing cap 
accelerometers, signals mould by more appropriately picked up 
non-contact type displacement transducers directly overlooking tbe 
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'..mft. Setter filters and amplifiers would also be needed. Shock- 
loads Clin be inflicted on the rotor in a more controlled manner by 
urjiiiQ electromagnets. The results are highly sensitive to data. 
Hence lutamatic data logging in the computer and processing are 

1 !ijp e r 3 1 1 V e . 

The rig is modular in design and can be readily employed to 
st-udy characteristics of composite shafts. Investigations and 
analytical modelling can be reformatted to study instability and 
define it in terms of acceptable levels of amplitude disturbances 
at different rotor speeds. 

The technique formulated in the present study needs to be 
implemented for fluid film bearings. This would be essential in 
order to establish its significance and practicability to real 


life rotors. 



59 


REFERENCES 

til Newl, ^ B.L., ’Shaft Whipping', General Electric Review, 
Vol. 27, 1924, p. 169. 

Cn ? imbaH, A.L., ‘Internal Friction Theory of Shaft Whipping', 
General Electric Review, Vol. 17, 1924, p. 244. 

C23 limball, A.L., ‘Internal Friction as a Cause of Shaft 
Whirling', Vol. 16, 1925, p. 724. 

C43 Ehrich, F.F., ‘Shaft Whirl Induced by Rotor Internal 
Damping', Int- of Appl. Mech., June 1964, p. 279. 

C51 Enrich, F.F. , 'The Dynamic Stability of Rotor/Stator Radial 
Rubs in Rotating Machinery', Jnl. Engg. for Ind. Trans. 
ASME, Vol. 91, No. 4, 1969, p. 1025. 

C63 Begg, I.C., 'Friction Induced Rotor Whirl - A Study in 
Stability', Jnl. Engg. for Industries, Trans. ASME, 1974, p. 
451 . 

171 Bansal, P.N. and Kirk, R.6. , 'Stability and Damped Critical 
Speeds of Rotor-Bearing Systems', Jnl. of Engg. for Ind., 
Nov. 1975, p. 1325. 

C83 Zorzi, E.S. and Nelson, H.D., 'Finite Element Simulation of 
Rotor Bearing Systems with Internal Damping', Jnl. Engg. for 
Ind., Trans. ASME, 1977, p. 71. 

191 Glasgow, D.A. and Nelson, H.D., 'Stability Analysis of 
Rotor— Bear ing Systems Using Component Mode Synthesis , Trans. 
ASME, Vol. 102, April 1980. 

C103 Tondl Ales, 'Some Problems of Rotor Dynamics', Chapman and 
Hall, London , 1 965 . 

C113 Rao, J.S., Rotor Dynamics : Wiley Eastern Limited, 1983. 

E123 Rao, J.S., and Gupta, K. , 'Introductory Course on Theory and 
Practice of Mechanical Vibrations', Wiley Eastern Limited, 
1984. 

C133 Rao, J.S., Gupta, K. and Vyas, N.S. , 'Blade Damping 

Measurements in a Spin Rig with Nozzle Passing Excitation 
Simulated by Electromagnets', Proc. Shock. Vib. Bull., US 
Naval Res. Labs., Proc. 56, Pt. 2, 1986, p. 109. 



40 


[141 J.S. and Vyas, N.S., 'Resonant Stress Calculation of A 

Turbine Blade with Damping as A Function of Rotor Speed and 
Vibratory Amplitude', ASMS 89-61-27. 

C153 Rao, J.S., and Vyas, N.S., 'Transient Stress Response of A 
Turbine Blade under Non-Linear Damping Effects', ASME 

90“GT~269. 



41 


APPENDIX 1 


The detailed SKpressions for computational process are as 
fol low.'-; ; 


G. 


1 


C IB-A) “+43 




'■(B-A) i X, + - Y, + (A+B) Y 4- (B^-AB+2)>.e 


pi pi 


12 — BT 

+ ^-(B-A) - X, - - Y, - (A+B) Y + (A*--AB+2)>.e 3.cos(T) 
p 1 p 1 o 

2 1 —AT 

+ C r.- X, - (B-A) - Y, - (B -AB+2) Y + (A+B>.e 
pi pi o 


^ 1 'I? — RT 

+ C- X, - (B-A) - Y + (A‘'-AB+2) Y +(A+B)>.e 3.sin(T)|-X 

pi pi o 


']-> 


or Cr- 

G„ = 


t 


1 

C <B-A)^+43 


(A.l) 


■ fc- X, + (B-A) - Y + {B^-AB+2) Y - ( A+B ) J . e 
[pi pi o 


1 ““BX 

+ -C- X, - (B-A) - Y, + (A -AB+2 > Y - (A+B>.e 3.cos(T) 
pi pi o 

t ^ ^ —AT 

+ C(B-A) - X, + - Y^ + (A+B) Y + (B''-AB+2) > . e 
p 1 p 1 o 

£(B-A) - X, + i Y, + (A+B) Y - (A^-A+B+2) > . e“®^3 -sin (T) |-Y 

o 


pi pi 


or 


^2 S' 




rCe>-A)‘*'^efj 

^ 2 (B-A) 

^ C (B-A) ^+43^ 


[c£zj 


(A. 2) 


dG. 






(B-A) +4 


TB+TA-1) i X - — Y + (-TA - TB+1 ) Y„ 
p 1 p 1 o 


o —AT 1 — BT 

+ (-TB + TAB-2T-B)> .e + £ Y + (2A-B)J.e 3.cos(T) 

p o 


+ cf- — X, + (TB-TA+1) - Y, +(TB^ - TAB+B^+2T) Y 

pi pi ° 
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m 

Mi 


I’ < A ' h '• 

r ih A5‘ ♦^r 
L, 


off 


X 

{(™i + Y + <2B-An.e 
P o 


(A. 3) 


♦ jaB-TA-l) I X + ^ Y + (TA + TB-1) Y 

p 1 p 1 ( 

« -BT 

4 <-TA“- 4 TAB -- 2T - A)}.e 3. cos (T> 

Y 1 —AT 

4 [ ' i ~ <2B-A>> Y + 1> -e 

p o 


11 a. + (TB-TA-l) I- Y- + (-TA^ + TAB-2T-A) Y 
pi P o 

e~®^3.sin(T)l 


4 (~TA~TB + in 


J 


(A. 4) 




(B”A) 


C (B-A)‘*+43'^ 


2 , 
<B-A) +4 


[*] 


2T 

p "1 


K- — X, + t-TB+TA-1) - + (-TB^-TA.B-2T-B> Y 

pi C 


+ CTA+TB-in.e 1 <2A-B) Y^ + l>.e ®^3.cos<T) 


2T 


+ CCC-TB+TA-U p + C-TA-TB+1) Y 


-AT 


+ <-TB^ + TAB - 2T-B>>.e 

+ £- 1 X + Y - <2A-B>> -e~®^].sin(T)] 

p 1 o j 


(A. 5) 



43 


dG, 

W 


2<A-B) 


C (B-A)'^+43^ 


H 


(B-A)^+4 


Cci Y, + (2B-A) Y -l>.e 
pi o 


+ C- — X, + (TB - TA - 1) - Y^ 

pi pi 

+ (-TA^+TAB-2T-A> Y + <-TA-TB+l ) > . e“®^3 . cos (T ) 

o 


1 —AT 

+ CC- X, + Y + (2B-A)J.e 
pi o 


+ t{-TB+TA+l) - X, + - Y, + (-TA-TB+1) Y 
p 1 p 1 o 

o —RT 

+ (TA^-TAB + 2T+A)>.e 3.sin(T) 


(A. 6) 
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APPENDIX 2 

A. The details of equipment used during experimentation are 
follows : 

1. Acce lerometers 

Make : B & K, Type 4370 

Charge Sensitivity : 10,10/10,20 PC " 

Frequency Range : up to 25 KHz. 

2. Charge Amplifier 

Make : B ?< K, Type 2635 

5 

Charge input : lO" PC Max 
Max Output ; 8V (mA) Peak 

Frequency Range i ACC/(0.2 or 2 to 100 KHz), 

Velo (1 or 10 to 10 KHz) 
disp. (1 or 10 to 1 KHz) 

3. Conditioning Amplifier 

Make : B K type 2626 

Charge input : 10^ pc max. 

Max. output : lOV (10mA) peak. 

Frequency Range : 0.3 Hz to 100 Khz. 

4. FFT Analyser Sm 2100B 

Make ; I watsu, Japan 


as 


No. of Channels : Two 

Sensitivity i ±0.1 to ±200 V/Full scale. 
Frequency Range : DC to 100 KHz. 
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5. Vibration meter 

Make : B & K Type 2511 

Frequency Range : 0.3 Hz to 15 KHz 
Accuracy : ±5% absolute 

6. Phase Meter 

Make : B & K, Type 2976 

No. of channels : Two 

Rotation freq. Range ; 2 Hz to 1 Khz. 
Photoelectric tacho. probe : Type MM 012 

7. Band Pass Filter 

i) Tunable band pass filter 
Make : B St K, Type 1621 
Input : IV Rms or 0.2 V Rms 
Pass bandwidth : Switchable 3*/* or 23*/. 
Frequency range : 0.2 Hz to 20 Khz. 

ii) Solid State Band pass filter. 

Make : KRQHN-HITE, USA, Type 3700<R> 
Input : ± 7V peak 
Pass bandwidth : 

Frequency Range = 0.2Hz to 20 KHz. 

Storage Oscilloscope 

Make : Kikusui, Type DSS 6521, Japan 

Sensitivity : 5 mV - 10 V/Div. 

No. of Channels : Two 
Sweep time mag : 5 times. 


8 . 
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9. Memory Hi Recorder 

Make : HIOK.I , Type 8801, Japan 
No. of Channels : Four 
Recording method : Thermal 
Sensitivity : 10 mV/Div to 5 V/Div. 

10. Stroboscope 

Make : General Radio, Type ; Strobotai 
Range : 110 to 150,000 RPM 


1538-A. 



The details of the experimen tal riq is as below : 
Motor 

Make : Tiillu-50 Water pump motor, 50 W 

Speed Range : 7000 rpm 

Base plate 

Material : Mild Steel 

Dimensions : im x 0.145 m x; 0.012 m 

Shaft 

Material : Stainless steel 
Length (is) : 0-5m 
Diameter (Id) : 0.01 m 

Modulus of elasticity (Es) =2.1 x 10” N/M'*“ 
Density ip) = 7800 Kg/wT' 

Disc 

Mass of disc 1 = 0.515 Kg 
Mass of dice 2 = 0.755 Kg 

Coupl ing 

Make : Luvjoy 

Bearings 

Make : SKF 

Outer dia- s 30 mm 


Inner dia. : 10 mm 
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APPENDIX 3 


RIGID ROTOR BALANCING 

The balancing procedures based on the influence coefficients 
measurement is as follows : 

i> Two convenient planes L and R and two measurement planes 
'a' and 'fa' (bearing locations) are chosen (See. Fig. AF.l). 

Taking and as the initial readings of vibration 

levels measured with phase angle and respectively at a 

desired rotational speed. 

ii) In the second run, a trial mass T placed at convenient 

location in plane R and observations and Z 6^ at 

measuring planes 'a' and 'b' were taken. Similarly in the third 
run, placing the trial mass T^ in plane L, z and R.^ z 

were measured. 

iii) The influence coefficient were obtained as 


=tR “ ‘”2 - 


“bL = 

“aR = ‘4 - 

4>/^r 

4i. = '4 


R^)/Tl 


L^)/Tl 


Hence 


f > 


- 

” 




, = - 

"^R 



«R 

^ 4 . 


. 4r 





The correction mass are obtained as 


4 “aL 

«R ' z r r — r— 

“bR “at ■ “aR “bL 
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ig. AF-1 The Block Diagram of Instrumentation Used in 
Rigid Rotor Balancing 







Fig. AF-2 


Balancing Measurements for Single Mass Rotor 






Fig. AF-3 


Balancing Heasure 


ments for Two Mass Rotor 




""l «aR - 4 

= 3— z Z 1— 

“bR "aL “aR %L 

The measured values are represented as vi*' shown in 

fiQ. AF.2 and fig. AF.3 for single mass rot«r niass 

rotor respectively- The phase angle measut were done 

taking positive in the direction of rotation Canti negative against 
the direction of rotation), referring the ma^^e^1 trial mass as 0 
in each correction planes. 

The correction masses and the residual unbala<^t:e obtained are 
given below, 

1. For single mass rotor ; 

= 3.294 gms ^ -143.6"^ 

= 0.389 gms ^ -134.6° 

L„ = 1.2 X 10~^ ^ 204° 

R 

= 0 . 6 >: 10 “" ^ 200 ^^ 

R 

2. For two mass rotor z 

FL = 0.887 gms ^ -98.6° 

R 

= 1.971 gms ^ 134.5° 

L„ = 0.0089 2 175° 

K 

= 0.0004 2 158° 

The rotors have always some residual unbalance, however well 


they might be balanced and are considered within tolerable limit. 
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